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DETERMINATION  OF  THE  LIGHT  ABSORPTION  INDEX  OF  SEAWATER  FROM  THE  LIGHT  FIELD 

PARAMETERS  OF  AN  ISOTROPIC  SOURCE 


[Pelevin,  V.  N.  and  T.  M.  Prokudlna,  Opredeleniye  velichiny  pokazatelya  pogloshcheniya 
sveta  morskoy  vodoy  po  parametram  svetovogo  polya  izotropnogo  Istochnika,  in:  Optics 
of  the  Ocean  and  Atmosphere  (Optika  okeana  i atmosfery) , Institute  of  Oceanology  of 
the  USSR  Academy  of  Sciences,  "Nauka"  Publishing  House,  Leningrad,  1972,  pp.  148-157; 
Russian] 

A relationship  between  the  light  absorption  Index^K.of  a medium  and  certain  para-^^' 
meters  of  the  state  of  optical  equilibrium,  permitting  the  calculation  of  K from 
measurements  of  the  light  field  in  a homogeneous  medium  at  large  depths,  was  estab- 
lished,bjtA.  A.  Gershun.l  This  method  of  determination  of  < has  also  been  treated 
in  theoretical  and  experimental  studies  by  other  investigators.^*^  A generalization 
of  this  metlUDd  to  the  case  of  a vertically  Inhomogeneous  plane-parallel  medium  uni- 
formly lllumYnated  with  natural  light  at  the  boundary  was  proposed  in  Refs.  4 and  5. 

The  latter  method  is  used  in  applied  marine  studies  conducted  by  the  Institute  of 
Oceanology,  lit  permits  the  measurement  of  the  vertical  distribution  of  the  absorption 
index  k(z)  in  the  sea  in  daylight  from  the  surface  down  to  depths  of  150-200  m. 


The  memod  of  determination  of  K from  the  parameters  of  the  light  field  of  an 
artificially  directed  light  source  in  a scattering  medium,  proposed  in  Ref.  6,  has  been 
fully  substantiated  theoretically,  but  the  procedure  Involved  in  obtaining  the  ex- 
perimental/data necessary  for  calculating  K is  very  complex.  In  the  simplest  case 
of  an  axlswmmetrlc  light  scurce,  it  is  necessary  to  move  the  irradlance  meter  in  the 
three  dir^tions  around  the  selected  point  in  the  course  of  the  measurements.  In 
addition/  experimental  data  on  spherical  irradlance  must  be  obtained. 


Tl*  irradlance  of  areas  placed  along  the  axis  of  the  light  beam  is  measured  /] 
with  T^rge  experimental  errors,  due  to  the  technical  difficulty  of  providing  for 
a satisfactory  cosine  characteristic  of  the  irradlance  meter  in  the  region  of  angles 
close  to  n/2.  For  .this  reason,  it  is  our  view  that  the  latter  method  is  not  suffi- 
ciently reliable  and  promising  for  the  determination  of  K in  situ  in  the  sea,  where 
the  [above-mentioned  measurements  are  particularly  difficult  to  perform.  The  organl- 
zatljon  of  systematic  measurements  of  absorption  index  < in  the  sea  at  night  and  at 
lari*  depths  in  the  daytime  remains  a very  pressing  problem. 

T!)e  present  study  sul)stantlates  the  metliod  of  determination  of  < from  the  para- 
meters of  the  light  field  of  an  isotropic  source  in  a homogeneous  medium.  The  in- 
struments and  method  of  processing  of  the  experimental  data  are  described.  Results 
of  determination  of  K are  given  for  several  regions  of  the  World  Ocean. 


js  I ight-scatte>i ng  ai 


Let  an  isotropic  llgl'.t  source  be  placed  in  a homogeneous  I ight-scatte>i ng  and 
absorbing  medium.  Tlic  boundaries  of  the  medium  arc  sufficiently  distant  from  each 
other.  It  is  well  known  that  the  spatial  emission  (absorption)  density  of  luminous 

] 

energy  at  some  point  of  tlie  light  field  is  equal  to  the  divergence  of  the  vector  H. 


ft 

Numbers  in  tlie  right  margin  indicate  pagination  in  tlie  original  text. 
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In  the  case  of  an  absorbing  medium,  this  relation  Is 

dlv  H - -ich,  (1) 

where  h = / Bdco  Is  the  spatial  Irradlance.  Integrating  expression  (1)  over  a closed 
4ir 

volume  located  between  spheres  with  radii  r and  r + dr  (with  the  center  of  the  spheres 
coinciding  with  the  light  source)  and  using  Gauss'  theorem,  we  obtain  the  equality 


where  S Is  the  sum  of  the  surfaces  of  the  two  spheres,  and  V Is  the  volume  Included 
between  them. 

Considering  that  as  a result  of  the  homogeneity  of  the  medium  and  Isotropy  of 

the  light  source,  the  vector  H Is  parallel  to  dS,  and  that  Its  modulus  H Is  constant 
over  the  entire  sphere,  we  rewrite  Eq.  (2)  as  follows: 

-4Kr-A/(/-)4-4-(r-f  rfr)5|//(/-)  + rf//|^-4.-r></rxA(r).  (3) 

Removing  the  parentheses,  discarding  second-order  terms  and  dividing  Eq.  (3) 
by  Anr^dr,  we  obtain 

1 dH 

Introducing  the  notation  h dr  ‘ we  arrive  at  an  expression  for  <:  J 


The  transport  vector  modulus  H may  be  expressed  In  terms  of  the  Irradlance  of 
areas  oriented  In  the  directions  toward  the  source  E and  away  from  the  source  E': 

H = E - E'  . 

It  should  be  noted  that  N.  G.  Boldyrev  once  proposed  a slmllaf  formula  for  deter- 
mining K in  the  light  field  of  an  Isotropic  source 

rfK  2 
£<tr  r • 

A more  rigorous  conclusion  leads  to  formula  (4).  It  is  easy  to  show  that  the  cooffi- 
2 

dent  a (r)  in  equality  (4)  may  be  expressed  In  terms  of  the  attenuation  index 

n V 

ot^(r)  of  the  flux  $(r)  through  a sphere  of  radius  r surrounding  the  Isotropic  source. 
By  definition 


whence 


Denoting 


U,V/{r), 
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we  obtain 


The  coefficient 


(8) 


is  the  "mean  cosine"  of  the  luminance  solid  at  a given  point  of  the  light  field. 

Thus,  the  absorption  index  of  the  medium  may  be  calculated  from  data  on  the  dis- 
tribution of  the  light  field  luminance  of  an  isotropic  source.  When  r formula 

(8)  becomes  the  known  relation  for  the  state  of  optical  equilibrium^  /151 


where  o is  the  "deep"  light  attenuation  index. 

The  light  field  characteristics  of  an  isotropic  source  in  the  sea  were  measured 
from  the  ships  AKADEMIK  VAVILOV  and  DKITRIY  MENDELEYEV.  An  outboard  photometric 
bench  designed  by  A.  S.  Suslyayev  was  used.  The  measuring  procedure  consisted  in 
the  following. 


Fig.  1.  Distribution  of  luminance  B(9)  in  the 
light  field  of  an  isotropic  source  at 
various  distances  r. 

St.  347,  southern  part  of  the  Sargasso  Sea, 

- power  of  the  source. 


An  Isotropic  light  source  (with  a quasl-sphertcal  light  intensity  Indioatrix), 
made  from  an  experimental  xenon  flash  lamp,  was  placed  on  the  upper  carriage  of  the 
bench  at  a depth  of  30  m.  The  photoelectric  detector,  a luminance  meter  with  a 
viewing  angle  of  1 or  10®,  made  from  an  FEU- 30  PM  and  covered  with  a ZhS16+S3S22 
glass  light  filter,  was  placed  on  the  lower  movable  carriage.  The  spectral  hand  of 
the  detector  was  495±35  nm.  Tlie  hinged  system  of  the  bench  made  it  possible  to  scan 


with  the  detector  in  the  plane  passing  through  the  source,  within  ±150*  of  the  dlrec-  r 
tlon  toward  this  source.  Measurements  of  the  luminance  solids  were  made  at  distances  1 1 
of  10  to  200  m between  the  source  and  the  detector.  The  received  signal  was  fed  to  /If  I | 
an  oscillograph  screen;  the  working  range  of  the  detector  was  adjusted  by  remotely  | | 
changing:  (a)  the  supply  voltage  of  the  PM;  (b)  the  load  resistance;  (c)  the  diaphragm  |j  | 
in  front  of  the  photocathode,  specifying  the  angle  or  reception  of  the  brightness  meter,  f | 
The  measurements  were  made  in  the  Black  Sea  and  in  various  water  areas  of  the  Atlantic  i | 
and  Pacific  Oceans.  The  vertical  homogeneity  of  the  water  was  checked  with  an  FPR  /15f  f | 
submersible  transparency  meter.  J | 


Fig.  2.  Irradiance  — - and  flux  through  a sphere  of  radius  r vs  distance  r 

^0  ‘**0 
from  Isotropic  source. 

Stations:  1 - 397,  2 - 367,  3 - 347,  4 - 342,  5 - 100,  6 - 410;  = 4ttJq  - power 

of  source. 

As  an  example.  Figure  1 shows  the  results  of  measurement  of  the  distribution 
of  luminance  B(0,r)  in  the  light  field  of  an  isotropic  source  In  the  Sargasso  Sea. 

Figure  2 shows  values  of  irradiance  E(r)  and  flux  'l’(r)  through  a sphere  of 
radius  r,  calculated  from  measurements  of  luminance  B(9,r). 

The  table  gives  values  of  K calculated  from  formula  (8),  and  also,  for  comparison, 
the  results  of  measurements  using  the  light  field  of  daylight, 4 made  on  water  samples 
with  a "Volna"  laboratory  differential  photometer,  and  using  the  signal  of  laser  pulse 
backscattcring  (data  of  D.  M.  Bravo-Zhlvotovskly  et  al.,  obtained  during  the  fifth 
cruise  of  the  ship  DMITRIY  MENDELEYEV).  The  last  column  of  the  table  gives  values 
of  the  mean  cosine  7],  used  in  calculations  of  K from  formula  (8).  Good  general  agree- 
ment between  the  results  of  measurements  using  different  methods  is  apparent.  A certain 
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Region  of  measurements 


Sargasso  Sea 

Central  region  of  Atlantic 
Ocean  (st.  342) 

North  tradewlnd  current  In 
the  Pacific  Ocean  (st.  367) 

Region  of  high  transparency 
near  Kuka  Islands  (st.  397) 

Pacific  Ocean,  region 
southeast  of  Japan  (st.  410) 

Black  Sea  (st.  100) 


Absorption  Index  K,  (values  given  with  decimal 
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discrepancy  of  the  data  may  be  due  not  only  to  Instrumental  errors,  but  also  to  dif- 
ferent conditions  of  measurement:  about  half  a day  elapsed  between  the  nocturnal  and 
diurnal  measurements,  and  the  ship  could  have  drifted  into  waters  with  somewhat  dif- 
ferent characteristics.  Measurements  of  the  absorption  Index  from  the  backscatterlng 
signal  give  the  average  value  of  K in  the  water  layer  from  the  surface  down  to  30-40  /15 
m,  whereas  the  light  field  of  the  Isotropic  source  is  measured  In  the  depth  range 
from  30  m and  below.  During  the  experiment,  certain  changes  are  possible  in  the  op- 
tical properties  of  the  sample  in  comparison  with  the  in  situ  state  - "spoilage  of 
the  sample"  - etc.  Thus,  the  results  of  the  comparison  of  tc  values  obtained  by  the 
proposed  method  with  the  data  of  other  measurements  should  in  our  view  be  considered 
satisfactory. 


Since  working  formula  (8)  was  derived  with  the  only  assumption  that  the  medium 
is  homogeneous,  the  error  in  the  determination  of  K (averaged  over  the  depth)  is  due 
solely  to  the  experimental  error  of  the  measurements  of  luminance  B(0) , which  were 
used  as  the  basis  of  the  calculation  of  the  absorption  index,  and  to  errors  of  the 
numerical  integration.  H,  h,  and  n wore  calculated  with  a "Minsk-22"  computer,  so 
that  the  integration  step  chosen  was  sufficiently  small,  and  the  calculation  error 
turned  out  to  be  much  (smaller  than  the  measuring  error. 

Obviously,  the  maximum  possible  relative  error  of  the  measurement  of  B(0)  is 


where  AB  and  AG  are  the  absolute  errors  of  measurement  of  luminance  and  of  the  angle 
to  which  this  luminance  Is  attributed. 


ab 

The  value  of  — Is  determined  by  the  operational  stability  of  the  source  and 

o 

by  the  calibration  accuracy  of  the  receiver,  and  In  our  case  amounts  to  5%.  The 
accuracy  of  adjustment  of  angle  A0  is  characterized  by  the  certificate  for  the  ro- 
tating devices  of  the  outboard  bench:  A0  = ±1*. 

9B 

Using  the  data  of  Fig.  1 for  determining  the  value  of  the  derivative  , we 

5B 

see  that  the  limiting  relative  error  of  the  luminance  measurement  — may  be  estimated 

B 

at  18X  in  the  range  0*  ^ 0 « 10*;  12%  for  10*  < 0 20®;  9%  for  20*  < 0 « 30*;  7% 

for  30*  < 0 < 60*;  6%  for  60*  < 0 ^ 180*.  On  the  basis  of  these  estimates,  and  con- 
sidering the  contribution  of  the  corresponding  sections  of  0 to  the  value  of  H,  we 

AH 

obtain  a limiting  error  ~ not  in  excess  of  10%.  The  limiting  error  for  was 

estimated  from  the  following  formula,  obtained  from  the  expression  for  the  total  dif- 
ferential a^: 

9 


• 


(10) 


where  rj^  and  T2  are  the  distances  between  which  the  mean  attenuation  index  of  flux  /155 
$ is  determined;  Ar  is  the  error  in  the  distance  measurements. 


Substituting  into  expression  (10)  the  mean  values  corresponding  to  our  measure- 
ments in  the  Black  Sea:  r,  = 50  m;  r„  = 150  m;  ot.  = 0.34  m“l;  Ar  = 0.25  m,  we  ob- 

12  9 

tain  a relative  error  of  determination  of  ct^  not  in  excess  of  7%.  For  measurements 
in  the  Sargasso  Sea,  rj^  = 50  m;  T2  = 200  m;  = 0.016  m and  as  a result. 


•10%  = 9%. 


We  will  calculate  the  limiting  error  in  the  determination  of  the  mean  cosine 
of  the  luminance  solid  7): 


^ a (0)  cosOsin®  rf#  ^ 

^ J D (e)sina</»  * 

4« 


(11) 


We  note  that  the  errors  in  B(0)  in  the  numerator  and  denominator  have  the  same 
sign  and  that  in  the  integration,  H and  h are  calculated  as  sums  of  the  type 
11  = la^B^;  h = where  a^  and  b^  are  coefficients  dependent  only  on  0.  We 

thus  obtain 
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Thus,  ^ — g — y is  the  weiglit  v/ith  which  the  error  in  enters  Into  the 

error  of  determination  of  the  mean  cosine  p.  Dividing  the  range  of  0 into  sections, 

°Bi 


as  indicated  above,  and  specifying  the  corresponding 


one  can  easily  calculate 


the  relative  error  in  the  value  of  the  mean  cosine  (12) . We  calculated  the  values 


6 


of  the  expression  In  parentheses  In  Eq.  (12)  from  data  of  specific  measurements  of 
luminance  solids  at  distances  of  20  and  100  m from  the  source  at  station  100  (Black 

Sea).  The  corresponding  values  of  •inoT  *-«  Ko  n anA  n nv 


lOOZ  were  found  to  be  0.65  and  0.8Z.  The 


calculation  also  showed  that  for  other  stations  as  well,  — ^ does  not  exceed  0.8%. 

n 

Aoc^  An 

From  expression  (8),  assuming  that  and  ~ are  Independent,  there  follows  /156 
an  expression  for  the  limiting  relative  error  of  determination  of  <: 

An  Aic 

Considering  the  above  estimates  of  and  -zr,  we  obtain  — 8%  in  the  Black 

, TT 

Ak 

Sea  and  — = <10%  in  the  Sargasso  Sea.  On  the  whole,  we  estimate  the  maximum  relative 
error  In  the  determination  of  the  absorption  index  at  not  more  than  10%. 


Actually,  the  total  error  Is  always  an  algebraic  sum  of  a large  number  of  par- 
ticular errors,  which  can  randomly  take  on  both  positive  and  negative  values  that 
partly  compensate  each  other.  Consequently,  it  may  be  expected  that  the  mean  statis- 
tical error  of  the  result  will  be  much  smaller  than  the  limiting  value  of  10%  given 
here. 


CONCLUSIONS 

1.  A method  of  determination  of  the  absorption  index  of  light  in  homogeneous 
waters  from  the  parameters  of  the  light  field  of  an  Isotropic  radiator  has  been  theore- 
tically substantiated  and  experimentally  verified. 

2.  In  situ  measurements  of  the  mean  absorption  index  in  the  495±35  nm  spectral 
range  were  carried  out  in  several  types  of  ocean  waters.  The  highest  measured  value 
of  K in  this  spectral  range  was  0.03±0,003  m“^  in  waters  of  the  Black  Sea;  the  loxjest 
value  v/as  0.014±0.001  m“l  in  the  high  transparency  region  near  the  Kuka  Islands  in 
the  Pacific  Ocean. 
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